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Abstract 

We present a full analytical solution for the localisation length in 
the one-dimensional Anderson model with weak diagonal disorder in 
the vicinity of the band centre. The results are obtained with the 
Hamiltonian map approach that turns out to be more effective than 
other known methods. The analytical expressions are supported by 
numerical data. We also discuss the implications of our results for the 
single-parameter scaling hypothesis. 
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1 Introduction 



Although it was introduced more than fifty years ago [Ij, the tight-binding 
model named after Anderson is still widely studied because it combines an 
elementary mathematical structure with non-trivial physical features. The 
simplicity of the definition notwithstanding, a complete analytical under- 
standing of the model is quite difficult to obtain, even in the one-dimensional 
(ID) case which is the most amenable to analytical treatment. Some funda- 
mental properties of the ID Anderson model, however, have long been known. 
In particular, it was rigorously proved that all eigenstates are localised in the 
ID Anderson model [2], unless the random potential exhibits specific spatial 
correlations (see [5] for a comprehensive treatment of localisation in models 
with correlated disorder). 

The localisation length is the key physical parameter which determines 
the spatial extension of the electronic states. A general formula for the local- 
isation length in the ID Anderson model is still not known, but expressions 
for the limit cases of strong and weak disorder have been derived. The weak- 
disorder case, in particular, can be studied using a perturbative approach 
due to Thouless [3]. Thouless' method gives a formula for the inverse local- 
isation length which works very well for most energies inside the band of the 
disorder-free model. This perturbative formula, however, has its fiaws: in 
fact, shortly after the publication of Ref. numerical calculations showed 
that Thouless' expression does not reproduce the correct value of the loc- 
alisation length at the band centre, i.e., for = [5]. A short time later, 
Kappus and Wegner were able to ascribe this discrepancy to a resonance 
effect which leads to a breakdown of the standard perturbation theory [B]. 
Almost at the same time, and without any explicit connection to the work of 
Ref. O |6] , an analytical expression of the localisation length for E = was 
derived in |7]. 

A thorough study of the "anomaly" at the band centre was eventually 
performed by Derrida and Gardner [8]. These authors showed that a "naive" 
perturbative approach was bound to fail not only at the band centre, but also 
for every other "rational" value of the energy, i.e., for E = 2cos(7rr) with 
r a rational number. Such an approach, in fact, gives an expansion of the 
localisation length which breaks down because some coefficients diverge. To 
avoid this pitfall, Derrida and Gardner devised a specific perturbative tech- 
nique which allowed them to analyse the anomalous behaviour of the localisa- 
tion length in the neighbourhood of the band centre. In particular, Derrida 
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and Gardner showed that the discrepancy between the numerical value of 
the localisation length and the prediction of Thouless' formula aroused be- 
cause the leading term in the correct expansion of the localisation length 
did not coincide with that obtained in Thouless' expansion. Derrida and 
Gardner also proved that a similar anomaly existed in the neighbourhood of 
E = ±2cos(7r/3) = ±1. This anomaly, however, manifested itself only in 
the next-to-leading term of Thouless' expansion and was therefore undetect- 
able within the second-order approximation. Derrida and Gardner surmised 
the existence of anomalies of a similar kind for all energies E = 2 cos(7rr) 
with r a rational number. This conjecture was later confirmed by the au- 
thors of Ref . , who focused their attention on the mathematical questions 
left open by Derrida and Gardner. Subsequently, in Ref. [101 |TT] a quasi- 
degenerate perturbation theory was developed with the aim of obtaining a 
uniform asymptotic expansion, in powers of the strength of the disorder, of 
the probability distribution for the ratio of the wavefunction amplitudes at 
neighbouring sites. Using this non-standard perturbative approach, the au- 
thors of Ref. [ini [H] re-obtained the results of Derrida and Gardner for the 
anomalies of the ID Anderson model at the band centre and the band edge. 

In the last few years the band-centre anomaly of the ID Anderson model 
has been considered from a new point of view, i.e., as a test case for the 
validity of the single-parameter scaling (SPS) hypothesis (see, e.g., [121 [131 
[11] and references therein). According to the SPS theory [TJl [IS], which 
is a cornerstone of our present understanding of Anderson localisation, the 
probability distribution of the conductance should depend only on a single 
free parameter. Studies of the ID Anderson model, however, have shown 
that this condition does not hold when the energy lies close to the band 
centre [121 [IS] ■ This has led to a renewed interest for the anomalies of the 
ID Anderson model in the neighbourhood of the band centre [121 [131 [S] ■ 

The purpose of this paper is to provide a complete analysis of the localisa- 
tion length in the neighbourhood of the band centre by means of the Hamilto- 
nian map approach. This formalism, introduced in [H], has proved to 
be quite an effective tool for the study of ID and quasi-lD disordered mod- 
els (see, e.g., [3j and references therein). The Hamiltonian map approach is 
based on the mathematical correspondence between the ID Anderson model 
and a classical parametric oscillator. This analogy makes possible to associ- 
ate quantum states of the Anderson model to phase-space trajectories of the 
parametric oscillator. In this way the phenomenon of quantum localisation 
can be understood in dynamical terms as energetic instability of a stochastic 
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oscillator Our first goal, therefore, is to use the Hamiltonian map ap- 
proach to derive in a mathematically simple and physically transparent way 
the results which have been obtained by Derrida and Gardner and subsequent 
authors with non-standard and intricate perturbative techniques. 

The second goal of this paper is to provide a detailed analysis of the 
transition from the anomalous behaviour of the model at the band centre to 
the normal regime away from E = 0. We focus our study on two physical 
magnitudes: the localisation length itself and the invariant distribution of 
the angle variable of the parametric oscillator (whose dynamics we describe 
in terms of action-angle variables). We find that both magnitudes exhibit 
a gradual crossover from the anomalous expressions they have at the band 
centre to the regular forms which are found for higher values of the energy. 
The slow progression of the invariant distribution towards a fiat form is im- 
portant because it calls into question the validity of the single-parameter 
scaling hypothesis. In fact, this hypothesis rests on the so-called random 
phase approximation, which is equivalent to the assumption that the invari- 
ant distribution of the angle variable be uniform [20] • As for the localisation 
length, our results confirm the blurred character of the transition which was 
recently pointed out in pi]. The analysis of Ref. [21], however, was based on 
purely numerical methods. Here we follow a different approach: we first ob- 
tain an analytical formula for the localisation length in the neighbourhood of 
the band centre and we then check its validity within the band by comparing 
its predictions with numerical results. 

The paper is organised as follows. In Sec. [2] we define the model and we 
describe how the Hamiltonian map approach works in the standard case. In 
Sec. |3]we show how the random- phase approximation fails in the neighbour- 
hood of the band centre and we derive the invariant distribution for the angle 
variable of the Hamiltonian map. In Sec. H] we use this result to obtain a 
general formula for the inverse localisation length. We draw our conclusions 
in Sec. |5l 
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2 Localisation length for non-resonant values 
of the energy 

We consider the ID Anderson model with weak disorder, defined by the 
Schrodinger equation 

tpn+l + tpn-l + £n4'n = Elp^- (l) 

Disorder is introduced in the model via the site energies En which are inde- 
pendent identically distributed random variables with zero average {Sn) = 
and variance 

(el) = cr' « 1. (2) 

Throughout this paper we will restrict our considerations to the weak-disorder 
case defined by condition ([2]). In this case the first two moments of the ran- 
dom site energies e„ provide sufficient information on the statistical prop- 
erties of the model ([T]). Our results, therefore, are vahd regardless of the 
specific form of the distribution of the site energies. 

The electronic states of the Anderson model ([1]) can be analysed in terms 
of the trajectories of a classical oscillator with Hamiltonian 

H = ^ + Imco'x' [1 + m] (3) 

where ^{t) is a frequency noise constituted by a succession of delta kicks of 
random strength 
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m= 2^ iJ{t-nT). 

n=—oo 

By integrating the dynamical equation of the kicked oscillator over the time 
interval between two kicks, one obtains the Hamiltonian map 

Xn+i = [cos {uT) - u^ri sin {ujT)]xn + —sin {ujT)pn /.^ 

muj [^) 
Pn+1 = —muj [sin (wT) + cos (wT)] x„ + cos {ujT) 

Here we have used the symbols Xn = x{t~) and pn = pit~) to represent the 
position and the momentum of the oscillator immediately before the n-th 
kick. After eliminating the momenta from the map one obtains the 
equation 

Xn+i + Xn-1 + ojS,n sin (coT) = 2 COS (ujT) Xn (5) 
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which has the same form of the Schrodinger equation ([T]) for the Anderson 
model provided that the following identities hold 



E = 2 cos {coT) (6) 

and 

£n = i^in sin {wT) . (7) 

Condition (|6]) corresponds to the dispersion relation for the Anderson model, 
if one makes the identification uT = /io linking the frequency u of the kicked 
oscillator and the period T of the kicks with the wavenumber fi of the elec- 
tronic states and the lattice step a of the model ([T]). Eq. ([7]), on the other 
hand, allows one to link the spatial disorder in the Anderson model ([T]) with 
the frequency noise of the parametric oscillator (j3]). 

The dynamical analysis of the kicked oscillator (|3]) becomes simpler if the 
Hamiltonian map (jl]) is written in action-angle variables introduced through 
the canonical transformation 



Xn. = \l sin 9n 

Pn = ^/2mU Jn COS 9n- 



One thus obtains 



sin 9n+i = [sin (6'„ + uT) - uj^n sin 9n sin (wT)] 
COS 9n+i = -=— [cos {9n + ujT) — uC,n siu 6'„ COS (ujT)] 



(8) 



with 



= ^ = 1 - 2wen sin 9^ COS e„ + sin2 9^. (9) 

Making use of the weak-disorder condition (|2]), the map expressed by Eqs. ([8]) 
and ([9]) can be written in the more compact form 



Qn+l = 9n + UjT + U^n siu^ 9^ + {u^nY slu^ 9^ COS 9n ,^ „x 

+ 0(^3) (mod27r) ^""^^ 
Jn+i = Jn{i - 2u;^„ sin 6'„ cos 9^ + u'^^l sin^ 6'„) . (11) 

In Eq. ffTOj) we have used the Landau symbol 0{a^) to denote neglected terms 
which, in the limit a — )■ 0, vanish like = {{eDY^"^ or faster (see, e.g., [22])- 
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The inverse localisation length (or Lyapunov exponent) is defined as 



lim 



1 ^ 



n=l 



ri-l 



N 



lim — lo£ 



n=l 



X'r 



In terms of the action-angle variables one can write the previous expression 
in the form 



1 ^ 1 
lim ^^logZ}„_i+ lim -log 



7V-s>oo 



n=l 



sin 6 



N 



sin 9() 



Except that at the band edge (where the angular variable tends to assume 
the values and vr), the second term in the right-hand side (rhs) of this 
identity vanishes; one is therefore left with 



1 ^ 



[12) 



n=l 



For weak and uncorrelated disorder one can expand the logarithm in the rhs 
of Eq. (JT2I) and factorise the disorder- angle correlators. Within the second- 
order approximation one thus obtains 



A = ^{ej [1 - 2(cos (2^„)) + (cos (4ft„))] + O {a') . 



(13) 



To proceed further, one must evaluate the averages of the trigonometric 
functions in the rhs of Eq. f lT3|) : it is therefore necessary to determine the 
invariant distribution p{6) for the angular map ffTOj) . 

Comparing Eqs. fflU]) and f lTT]) . it is easy to see that the angle variable 
has a faster dynamics than the action variable and that a limited number of 
map iterations suffices to make the 6 variable sweep the whole [0, 2n] interval. 
One can therefore expect that the distribution for the angle variable should 
quickly reach a uniform invariant form, 



(14) 



We remark that the assumption that the invariant distribution has the flat 
form corresponds to the so-called random phase approximation [T2| 120]. 
This approximation assumes that, for a random sample of length L, the 
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phases of the reflection and transmission amphtudes are uniformly distrib- 
uted over [0, 27r] when the length of the sample is much larger than the 
localisation length, L ^ 1/A. The random phase approximation is usually 
invoked in the SPS theory of localisation [121 H]; when condition f lT^ is not 
met, therefore, one may expect the SPS theory to fail. 

Besides being a key ingredient in the derivation of the SPS theory, con- 
dition f[T^ is also required to obtain Thouless' formula. In fact, Thouless' 
expression for the localisation length can be easily recovered from Eq. (1T3|) 
by computing the averages in the rhs of this equation with the invariant 
distribution (fT^ . In this way one arrives at the formula 

The analysis of the map (fTU|) shows that the assumption (IT^ is usually 
correct; however, it fails whenever the parameter uT is a rational multiple 
of TT. For uT = vr/g, in fact, the angular map f fTOj) has periodic orbits of 
period 2q in the absence of disorder. When weak disorder is introduced, the 
resulting orbits linger around the unperturbed periodic orbits; as a result, 
the invariant distribution ceases to be fiat. The modulation of the invariant 
measure is the origin of the failure of Thouless' formula (fT5|) at the band edge, 
i.e., for uT = or uT = ±7r, and at the band centre, i.e., for uT = tt/2. 
In general terms, the invariant distribution acquires a modulation for all the 
"rational" values of the energy E = 2cos(7rp/g); however, within the limits 
of the second-order approximation, the localisation length differs from the 
standard value predicted by Eq. f|T5|) only for |g| < 2, i.e., at the band edge 
and at the band centre. Higher values of q, in fact, produce a high-frequency 
modulation of the invariant distribution which does not affect the averages 
of the low-order harmonics in the rhs of Eq. flT^ . 

The anomalous behaviour in the neighbourhood of the band edge was 
thoroughly analysed with the Hamiltonian map approach in [18]. We now 
turn our attention to the neighbourhood of the band centre, devoting the 
next section to the derivation of the invariant distribution in this region. 



3 The invariant distribution 

We consider the case in which 

u;T=^ + 6 (16) 
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with 5—7-0. Taking into account the dispersion relation 
Eq. (fTB]) imphes that 

E = -2sin5 ~ -25, 



for 5 — > 



i.e., the energy hes in a neighbourhood of the band centre. From Eq. f lTU]) it 
is easy to see that, when condition f lTB]) is met, the unperturbed angular map 
has almost-periodic orbits of period 4. This implies that 9n+4 must lie close 
to 6n even when disorder is present. We can compute the difference 6'„+4 — 9n 
by considering the fourth iterate of the angular map ( ITOj) . Neglecting terms 
of the form ^„5, which are of order 0{a6), and making use of the statistical 
independence of the variables one can write the fourth iterate of the 
angular map ffTOj) in the form 



n+4 



^„ + 45 + - + in+2) [1 



COS (20„)] 



I {in+1 + a+s) [1 + COS (2^„)] - y sin (4^, 



(mod 2-7?) 



(17) 



The sign of approximate identity in Eq. f[T7j) is due to the fact that the rhs is 
a truncated expansion in both the a and 5 variables in which we omit O (cr^) 
terms of order higher than the second in the disorder strength, second-order 
terms O (5^) in the energy shift 5, and cross-terms 0((j5). Note that, in the 
limit case 5 = 0, Eq. (fT7|) reduces to the form used in [18] to analyse the 
exact band-centre case. By going to the continuum limit one can replace the 
map (1171) with the Ito stochastic differential equation 



d^ 



+ 



45 + y sin (4^) 



0" 




dt+ \ —[l-cos{29)]dWi 




y [1 + cos {26)] AW2 



where Wiit) and W2{t) are independent Wiener processes. The Ito stochastic 
equation ([T5|) . together with an initial condition Oito) = 9q, completely de- 
termines the stochastic process 9{t). 

As is well known [23], the 6{t) process can be specified also in terms of 
the conditional probability distribution 

p = p{9,t\9o,to) 

which is obtained by solving the associated Fokker-Planck equation 



dp 
di 



a 
T 



d_ 



1647 + 2 sin (4^) 1 P 



+ 



86' 



[(3 + cos(4^))p] ^ (19) 
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with the initial condition p {9, to\9o, to) = 5{6 — 9q). 
It is convenient to introduce the parameter 



25 E 

"=-1^^^ (20) 

which measures the distance of the energy from the band centre on a scale 
defined by the strength of the disorder. From a mathematical point of view, 
the parameter (120]) represents a ratio which is held fixed while the double 
limit 

E ^0 and ^ 

is taken. Inserting the new parameter in Eq. (fT9|) and rescaling time accord- 
ing to the identity 

one can cast the Fokker-Planck equation (fT9!) in the form 

1^ = ^ { [8x + 2 sin(4^^)] P + ^ [(3 + cos(4^)) p] } . (21) 

The stationary solution of the Fokker-Planck Eq. fl2Tl) is the invariant 
distribution p{9) we are interested into. Therefore we must find the solution 
of the stationary equation 

^ I [8x + 2 sin m P (^) + ^ [(3 + cos (4^)) p (9)] | = (22) 

which satisfies the conditions of periodicity 

p{9 + 27r)=p (9) (23) 



and normalisation 



2tt 

p{9) d9 = 1. 



Integrating once Eq. (1221) one obtains the first-order ordinary differential 
equation 

dp _2 sin (4^) - 8h C 

d9 ~ 3 + cos (4^) ^ ^ 3 + cos (4^) 
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with C being an integration constant. The general solution of this equation 
has the form 



P{0) 



with 



m 



2Q-s>cF{e) 
V^3 + cos(4^) 



d0 



P(0) + 



C 



,8xF{0) 



v^3 + cos (40) 



d(j) 



(24) 



3 + cos (40) 



I 



4^2 I 



Tcn + arctan 



tan{29) \ for n-<9<{n+ 1)-. 
v2 J J 2 2 

(25) 

In Eq. ( l25l) n = 0, ±1, ±2, . . . and arctan denotes the inverse of the tangent 
function having the interval [0, vr] as domain. The solution contains the 
constants C and p(0) which must be determined from the periodicity and 
normalisation conditions. The periodicity condition ( 123|) implies that 



C = 2p(0) (^e^^'^" ~ // 



2n g8>fF(</.) 



v/3 + cos(40) 



d0 



Substituting this identity in Eq. ( l24l) and imposing the normalisation condi- 
tion, one eventually obtains that the invariant distribution has the form 



piO) 



1 



Ni>c) ^3 + cos(4^) 



v^3 + cos(40) 



where the normalisation factor N{k) is equal to 



+ 



27r 



,8>fF(0) 



e v/3 + cos(40) 
(26) 



N{x) 



2-n 



^/3Tms{Ae) 



^4^271 



g8^F(<A) 



+ 



27r 



,8^F(</.) 



/o v^3 + cos(40) Je y3Tcos(40) 

(27) 

To the best of our knowledge, the general expression (12^ of the invariant dis- 
tribution has not been discussed before. Ref. [llj contains equations which, 
if combined, are equivalent to Eq. f l26|) but the authors of Ref. |TT] worked 
out the explicit form of the invariant distribution only for the limit cases of 
|x| ^ 1 and |x| ^ 1. The limit forms of the invariant distribution are much 
easier to use (see below), but in this paper we prefer to focus our attention 
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on the general expression (12^ . because it is the only distribution that al- 
lows one to analyse the behaviour of the localisation length over the whole 
neighbourhood of the band centre. 

An important property of the invariant distribution fl26|) is that it is 7r/2- 
periodic, i.e., 

p{0 + l)=pi0). (28) 

In physical terms the periodicity condition ( 125]) should be seen as an expected 
manifestation of the fact that unperturbed orbits of the map fllUp have period 
4 for E = 0. From a mathematical point of view, Eq. f l28|) can be obtained 
by first noting that the definition fl25|) implies that 



F 



F{9) + 



vr 



4^2' 



(29) 



Keeping in mind Eqs. (!25|) and (129|) . one can analyse the behaviour of the 
integral terms in the representation of the invariant distribution. It is 
easy to show that 



/ 

Jo 



/3 + cos(40) 



7r/2 



+ 



^3 + COS (40) 



(30) 



v^3 + cos(40) 



and 



2-IT 



,?,kF{4,) 



e+TT/2 a/3 + cos(40) 



,4^271 



2tt 



o8^F{<f,) 



e J 3 + cos (40) 



d0 



tt/2 



(31) 



v^3Tcos(40) 



d0. 



Making use of the identities ( 129|) . (!30|) . and (I3T]) . one can evaluate p{9 + Tx /2) 
from Eq. (!26ll and thus obtain the periodicity condition (l28l) . 

The validity of the analytical expression (l26ll is corroborated by numerical 
computations, as can be seen from Figs. [1] and [21 In Fig.[T]the theoretical dis- 
tribution fl26l) is plotted together with the distributions which were obtained 
numerically for several values of x. As can be clearly seen, the invariant 
distribution ( l26l) matches very well the numerical data. Fig. [T] confirms that 
the invariant distribution is not uniform in the neighbourhood of the band 
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Figure 1: (Colour online) Invariant distribution p{6) for various values of x. 
The points represent the numerically computed invariant distribution, while 
the lines correspond to Eq. fl26|) . 



centre. The amplitude of the modulation is largest at the exact band centre, 
i.e., for X = 0, and gradually decreases as as the energy moves away from 
the band centre (i.e., for increasing values of 

The numerical data displayed in Fig.[T]were obtained for disorder strength 
= 0.01, but the results are the same for different values of cr^. This is 
made evident by Fig. [2], which represents the invariant distribution of the 
map ffTOj) at the band centre {E = 0) computed for different values of the 
disorder strength. The good data collapse on the line corresponding to the 
theoretical form ( l26l) confirms another feature of Eq. ( l26i) . i.e., that in the 
neighbourhood of the band centre the invariant distribution is determined 
by the parameter x and does not depend on the disorder strength. The 
data in Fig. [2] were obtained for x = but the same collapse occurs also for 
non- vanishing values of x. 

As mentioned before, the general expression (1261) for the invariant dis- 
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tribution reduces to simpler forms in the limit cases |x| -C 1 and |x| ^ 1, 
which are the only ones considered in the literature [H [11] . We devote the 
rest of this Section to the analysis of these particular cases. 



3.1 The invariant distribution for |>z:| <C 1 

In the limit case |x| ^ 1 one can expand both the expression fl26l) and the 
normalisation factor ( 1271) in powers of h. Carrying out the calculations, one 
eventually obtains that the limit form of the invariant distribution is 



2p(0) 



^3 + cos(4^) 



X 



£(l) 
r(i) 



^3 + cos(40) 



d0 



+ 0(> 
(32) 
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with the constant factor 2p(0) given by 



2p(0) 




+ 0( 



Eq. ( !32l) coincides with the corresponding result first obtained by Derrida and 
Gardner [8] ; the apparent difference is due to the election of these authors of 
a complex representation for the function ( 125|) . 

At the band centre, i.e., for x = 0, expression (!32|) reduces to the form 



2r(!) 



^r(i) v/3 + cos(4^) 2K (^) v/3 + cos(4^)' 



In the previous equation the symbol K(A;) stands for the complete elliptic 
integral of the first kind. 



K(A;) 

and we used the identity 



n/2 



1 — fc^ sin^ I 



8r(|)- 



3.2 The invariant distribution for |>f| 1 

We now consider the case |x| ^ 1. In this limit, the asymptotic expansion 
of p{6) is obtained with less effort taking as a starting point the differential 
equation ( 122|) . rather than the general form (126|) of the invariant distribution. 
After writing Eq. (!22|) as 



p{d) - sin(4e)p(^^) + ^ (3 + cos(4^)) % 
4x ox au 



0, 



one can expand the solution in powers of 1/x, 

X X- 



(33) 



(34) 
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Substituting the expansion 
equations 



in Eq. fl33p one obtains the hierarchy of 



d9 



and 



d 



P 



(n) _ 



sin (461) p 



(n-l) 







1 r/rt^""^) 
-(3 + cos(4.))Ji^ 







(35) 



(36) 



with n = 1,2,... For the expansion to satisfy the condition of period- 
icity and normahsation, the solution of Eq. fl35|) must be 27r-periodic and 
normahsed, while the solutions of the higher-order equations (1361) must be 
27r-periodic and satisfy the conditions 



2tt 



p(")(^)d^ = 0. 



In this way one can easily obtain that the behaviour of the invariant distri- 
bution for X 1 is 



PW = ^ + -^sin(4^^)-AT|- 
Ztt xStt x"^ 167r 



cos(40) 



cos(86 



+ 0(^). (37) 



Eq. (13 7p coincides with the result previously obtained in [S]. Note that as |x| 
increases, i.e., as the energy moves away from the band centre, the invariant 
distribution tends to recover its flat form. The uniform limit, however, is 
reached with power-law convergence: this implies that the transition from 
anomalous to regular behaviour is smeared out over a wide energy range. 



4 The localisation length 

Having obtained the invariant distribution f l26|) . one can compute the rhs 
of expression f[T^ . After observing that the average of cos(26') vanishes 
because the invariant distribution has period 7r/2, one obtains that the inverse 
localisation length is equal to 



A 



X 



8(l + f) 



1 + 



2n 



cos(40) 



iv(x) Jo VsT^J^^^i^ 



-.d(p 



a/3 + cos(40) Je a/3 + cos(40) 



2n 



-.d(f) 



(38) 



d^ 
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The general expression ( 155]) is the central result of this paper. Its absence 
from previous works on the band-centre anomaly is probably explained by 
the fact that in the literature only the limit cases |x| -C 1 and |x| 3> 1 have 
received attention so far. We discuss these limit cases below; we would hke 
to emphasise, however, that only the general expression fl38|) makes possible 
to describe the behaviour of the localisation length over the whole range of 
the parameter x. 

A remark is in order here. Because the invariant distribution (!26|) was 
obtained under the assumption E ~ —26 — )■ 0, in the energy range where 
Eq. (138|) is rigorously valid, it is equivalent to the following form 



which is obtained by setting E = in Eq. (155]) . However, as discussed below, 
we have found that Eq. fl55]) has a validity of its own because it provides 
a very effective interpolation between formula ( 139]) . which is applicable in a 
neighbourhood of the band centre, and Thouless' expression (fT5]) which can 
be used in the rest of the band (with the exception of the band edges). 

It is not difficult to evaluate numerically Eq. f l39p : this operation makes 
possible to compare analytical and numerical results for the inverse localisa- 
tion length. As shown by Fig. |3l the theoretical predictions of Eq. f l39]) fit 
very well the numerical data in a neighbourhood of the band centre. For the 
case represented in Fig. [3] the disorder strength was set at = 0.01; the 
considered energy interval [—0.1, 0.1] corresponds to values of x in the range 
[—10, 10]. It is easy to notice that the actual localisation length exhibits a 
non-negligible difference from the values predicted by Thouless' formula even 
for energies lying quite away from the band centre. This is a consequence of 
the slow convergence of the invariant distribution (137]) to the uniform limit 
discussed in Sec. 13. 2j we can conclude that the phenomenon of anomalous 
localisation is not restricted to an infinitesimal neighbourhood of the band 
centre, but can be detected over a finite interval of energy values. 

A close examination of Fig. [3] shows that, for energies l^'l > 0.05, the nu- 
merical data tend to overlap with Thouless' formula, while Eq. f l39]) fiattens. 
This is due to the fact that, as x increases, Eq. fl39]) tends to a constant limit 
(see Eq. fllT]) below). The physical reason is that Eq. f l5I?]) is valid only as 



A 




(39) 
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Figure 3: (Colour online) Inverse localisation length A versus energy E. The 
dots represent numerically computed values of the Lyapunov exponent; the 
dashed line is given by Eq. fl39|) : the dotted line corresponds to Thouless' 
formula. 



\E\ — )■ 0. It turns out, however, that Eq. ( l38l) works very well over the whole 
energy band. This is born out by Fig. HJ which compares the numerical data 
both with Thouless' formula ( !T5|) and with Eq. ( !38|) . That Eq. ( 138|) works so 
well over the whole energy band is not surprising. In fact, the rhs of Eq. (138!) 
is the product of Thouless' formula and of a corrective factor. The latter 
is responsible for the anomalous behaviour of the localisation length in the 
neighbourhood of the band centre, but tends to unity for increasing values of 
In this way Thouless' expression is recovered away from the band centre 
and the anomaly is taken into account by the corrective factor. 



18 



0.0017 
0.0016 
0.0015 
0.0014 
0.0013 
0.0012 
0.0011 



0.00125 



\ ^ 0. 

\ 



00120 




\ 0.00115 
\ -0.2 



-1 



-0.5 



I 




/ 



Numerical data 
Analytical formula 
Thouless formula 





E 



0.5 



/ 



Figure 4: (Colour online) Inverse localisation length A versus energy E. The 
dots represent numerically computed values of the Lyapunov exponent; the 
dashed line is given by Eq. fl38|) : the dotted line corresponds to Thouless' 
formula. The inset shows in more detail the same data in the neighbourhood 
of the band centre. 



4.1 Localisation length for |x| <C 1 

We now turn our attention to the limit cases |x| ^ 1 and |x| ^ 1. In the 
case |x| -C 1, i.e., when the energy lies very close to the band centre, one can 
evaluate the rhs of Eq. (JT3ll with the approximate invariant distribution (l32ll . 
In this way one obtains 



£(!) 



[l + 0(x)] 



(40) 
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in agreement with the known result [HI [HI IB] • We remark that Eq. fl4UI) can 
also be written in the equivalent form 



A 




1 [l + 0(x)] 



which is sometimes used in the literature [3123]. The symbol E(fc) represents 
the complete elliptic integral of the second kind, i.e., 



4.2 Localisation length for |x| 1 

In the opposite case, i.e., for |x| ^ 1, the energy, although close to the band 
centre in absolute terms, moves away from it on the energy scale set by the 
strength of the disorder. When |x| ^ 1 the invariant distribution takes the 
form fl371) and the inverse localisation length fl39|) becomes 



Once more, the result coincides with the formula first derived by Derrida and 
Gardner p[II]. 

In the limit |x| — t- oo the inverse localisation length tends to the value 
predicted by Thouless' formula f|T5l) : as already remarked the transition to the 
regular limit is not sharp because of the power-law decay of the anomalous 
term in the rhs of Eq. f HTj) . It is important to remark, however, that the 
limit |x| — )■ cxD in Eq. f l41l) can be taken only by letting the intensity of the 
disorder cx^ tend to zero faster than the energy E. Simply increasing |x| 
while keeping fixed the disorder strength eventually leads to the breakdown 
of the formula (HTj) . which was derived under the assumption that E 0. 
This can be seen in Fig IH) for energies \E\ > 0.05 the inverse localisation 
length tends to move away from the anomalous expression ( l39l) and closer to 
Thouless' formula. 





(41) 
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5 Conclusions 



In this work we have analysed the anomalous localisation of the eigenstates 
of the ID Anderson model for energies close to the band centre. Using the 
Hamiltonian map approach, we have derived two main results: the general 
expression ( l26ll of the invariant distribution for the map ( JTOl) and the cor- 
responding formula f l38|) for the inverse localisation length. The first result 
shows that the random phase approximation, which is an essential ingredi- 
ent of the single-parameter scaling theory, fails close to the band centre. The 
failure is not restricted to a negligible energy range because the convergence 
of the distribution (^E^ to the uniform limit follows a power-law behaviour. 

The invariant distribution fl2^ also allowed us to obtain the analytical 
expression (!38|) of the inverse localisation length. We have analytically proved 
that the general formula (!38|) describes the Kappus-Wegner anomaly around 
the band centre; because it reduces to Thouless' formula when the energy 
moves away from the band centre, Eq. (1381) works well over the whole energy 
band, as confirmed by the numerical data. 

This paper provides new insight for a critical investigation of the SPS 
theory with the Hamiltonian map formalism. In this work we have focused 
our attention to the ID map ffTOl) for the angle variable and in this way we 
have been able to show that the random phase approximation, and hence 
the SPS theory, fail in a finite energy interval around the band centre. How- 
ever, an extension of our analysis to the complete action-angle map given by 
Eqs. (fTO!) and (fTT!) is required in order to study the statistical properties of 
the conductance of a random segment. In fact, the conductance is related to 
the action variable, which is statistically correlated with the angle variable: 
hence a simultaneous study of both variables cannot be avoided. We plan to 
address this issue in a future work. 
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